Energy density for chiral lattice fermions with chemical potential 
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We study a recently proposed formulation of overlap fermions at finite density. In particular we 
compute the energy density as a function of the chemical potential and the temperature. It is shown 
that overlap fermions with chemical potential approach the correct continuum behavior. 
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I. INTRODUCTION 

Over the last two decades lattice gauge theory was 
turned into a powerful qualitative tool for analyzing 
QCD. This progress is in part due to the advances in 
algorithms and computer technology, but also on the con- 
ceptual side important breakthroughs were made. Most 
prominent among these is the correct implementation of 
chiral symmetry on the lattice based on the Ginsparg- 
Wilson equation for the Dirac operator 

An application of lattice techniques which has seen a 
lot of attention in recent years, is the study of QCD at 
finite temperature. The lattice implementation of the 
chemical potential /i, necessary for such an analysis, is 
not straightforward, however. It is well known [2j, that 
a naive introduction leads to /i 2 /a 2 contributions which 
diverge in the continuum limit when the lattice spacing 
a is sent to zero. For more traditional formulations, such 
as the Wilson or staggered Dirac operators, the problem 
has been solved by introducing the chemical potential in 
the same way as the 4-component of the gauge field. 

A satisfactory implementation of the chemical poten- 
tial should be compatible with chiral symmetry on the 
lattice based on the Ginsparg- Wilson equation. When 
attempting to introduce the chemical potential into the 
only solution of the Ginsparg- Wilson equation know in 
closed form, the overlap operator a potential prob- 
lem quickly surfaces: defining the sign function of a non- 
hermitian matrix. In [|[ Bloch and Wettig proposed a 
solution based on an analytic continuation of the sign 
function into the complex plane. It was shown, that the 
eigenvalue spectra of this construction match the expec- 
tations from random matrix theory. 

In this letter we analyze the proposal [4J further 
and study the energy density of free, massless overlap 
fermions with chemical potential. The dependence of the 
energy density on fi and the temperature T allows for a 
detailed analysis of the lattice formulation at finite den- 
sity. Of particular interest will be the question whether 
the analytic continuation of the sign function produces 
divergent /J? /a 2 terms. Our study indicates the absence 
of such contributions and we find that the \i and T de- 
pendence of the energy density is approached correctly. 



II. SETUP OF THE CALCULATION 

The overlap Dirac operator D(fi) for fermions with a 
chemical potential fi is given as 
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D(n) = ~ [l-7 5 sign , 
a 

H(jjl) = 75 [1 — aDw(fJ-)] ■ 
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The sign function may be defined through the spectral 
theorem for matrices. Dw(^) denotes the usual Wilson 
Dirac operator, 
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For later use we distinguish between the lattice spacing 
a in spatial direction and the temporal lattice constant 
0,4. Periodic boundary conditions are used in the spatial 
directions, while in time direction we apply anti-periodic 
boundary conditions. The chemical potential // is cou- 
pled in the usual exponential form Q ■ 

For vanishing /1 the Wilson Dirac operator is 75- 
hermitian, i.e., 75-Diy (0)75 = D w(0) ■ This implies that 
H (0) is a hermitian matrix. As soon as the chemical po- 
tential fj, is turned on, 75-hermiticity no longer holds, and 
iJ(/x) is a non-hermitian, general matrix. This fact has 
two important consequences: Firstly, the eigenvalues of 
are no longer real and the sign function for a com- 
plex number has to be defined in the spectral representa- 
tion of sign H([i). Secondly, the spectral representation 
has to be formulated using left and right eigenvectors. 
This latter problem will be dealt with later when we dis- 
cuss the evaluation of sign H (fi). For the sign function 
of a complex number we use the analytic continuation 
proposed in Q and define the sign function through the 
sign of the real part 



sign [x + iy ) — sign (x) 



(3) 
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The observable we study here is the energy density 
defined as 
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Here 7i is the Hamiltonian of the system, M denotes the 
number operator and (3 = 1/T is the inverse temperature 
(in our units the Boltzmann constant k is set to k = 1). 
The derivatives in the second line are taken such that 
/3/i = c = const. 

The continuum result for the subtracted energy density 
of free massless fermions reads (see, e.g., @) 



e(ji) - e(0) 



4tt 2 



2 rji'2 



+ 2 



(5) 



When working on the lattice, the inverse tempera- 
ture P is given by the lattice extent in 4-direction, i.e., 
/3 = N 4 a 4 . Thus the derivative d/d/3 in ^ turns into 
N^ 1 d/da 4 . The partition function Z is given by the 
fcrmion determinant det D which we write as the prod- 
uct over all eigenvalues A n . We thus find 
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III. EVALUATION OF THE EIGENVALUES 

According to (O, for the evaluation of e(fi) the eigen- 
values A„ of the Dirac operator D have to be computed. 
This is done in three steps: First we bring the Dirac oper- 
ator for free fermions to 4 x 4 block-diagonal form, using 
Fourier transformation. Subsequently the spectral repre- 
sentation is applied to the 4x4 blocks of H to evaluate 
sign H . Finally the eigenvalues of the blocks of D are 
computed and by summing over the discrete momenta 
all eigenvalues are obtained. 

Following this strategy, one finds for the Fourier trans- 
form H of H, 



H = 75/15 + i'ys'yi'YvK 



(7) 



with 



h>5 = 1 — — cos(apj)] — — [1 — cos(a4(p4 — iy))} , 

■ , a 4 
3=1 

hj = — sin(apj) for j = 1,2,3, 

hi = sm{a 4 {p 4 — ifj,)) . (8) 

(24 



The spatial momenta are given by pj = 2irkj/aN, 
where N is the number of lattice points in the spatial 
directions and kj = 0, 1 ... N — 1. The momenta in time- 
direction are p 4 — n(2k 4 + l)/a 4 N 4 , k 4 — 0,1 ... N 4 — 1. 

The remaining diagonalization of H is similar to the 
construction of the left- and right-eigenfunctions for the 
free Dirac operator. One finds that H has two different, 
doubly degenerate eigenvalues 



■ s , a<3 = a 4 



\jw + h 2 , (9) 



where h 2 = ^2 v h 2 . The corresponding left- and right- 
eigenvectors, lj and Tj are given by 

h = lf ] [H + sl] , l 2 = lf\H + st\ , 



ri = [H + sl]rj 0) , r 2 = [H + sljr^ 
r 3 = [H - sl]4 0) ,n = [H- sl]r { 4 0) 



The constant spinors ly , r,- are (T is transposition) 
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4 0) = rf r = c(0,0,l,0),4 0) =rf T = c(0,0,0,l). 

(11) 

The constant c = (2s(s + h^))" 1 / 2 ensures the correct 
normalization, such that the eigenvectors obey Urj = Sij. 

Using these eigenvectors and the spectral theorem we 
find for sign H the simple result 



sign if = ^r>ign(A ; ir,/, = *&Qh 



(12) 
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Plugging this back into the overlap formula ([T]) and di- 
agonalizing the remaining 4x4 problem one finds two 
different eigenvalues for the overlap operator at a given 
momentum, 



sig 



n(Vh 2 



h 2 )hr Q ± iVh 2 



s/h? + h 



(13) 



where each of the two eigenvalues is twofold degenerate. 
The momentum dependence enters through the compo- 
nents h v , /15 defined in (JHJ) . In the spectral sum (J6j> the la- 
bel n runs over all momenta and the eigenvalues at fixed 
momentum as given in (|13p . The necessary derivative 
with respect to a 4 is straightforward to compute in closed 
form, and the spectral sum ([B]) can then be summed nu- 
merically. The argument of the sign function cannot be- 
come purely imaginary on a finite lattice, and no <5-like 
terms occur. We remark, that after taking the derivative 
with respect to a 4 , we set a — a 4 — 1, i.e., all the results 
we present are in lattice units. 
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FIG. 1: The energy density e(fi) — e(0) as a function of /i 4 (all 
in lattice units). The symbols (connected to guide the eye) 
are for various lattice sizes, the dashed line is the continuum 
result. 



IV. RESULTS 

We begin the discussion of our results with Fig. 1, 
where we show the subtracted energy density e(/i) — e(0) 
as a function of /u 4 for three different lattice volumes. For 
those lattices all 4 sides have equal length, i.e., in the 
thermodynamic limit they correspond to zero tempera- 
ture. Thus, according to ([5]), we expect the data (sym- 
bols in Fig. 1) to approach the continuum form /i /Air 2 
(dashed line) as the 4-d volume is sent to infinity. 

The figure clearly shows that the lattice data are pre- 
dominantly linear when plotted versus ll and that for 
small \i they approach the continuum curve when the 
volume is increased. It is, however, obvious that also on 
our largest lattice still a discrepancy remains for larger 
ll. In particular one finds a slight curvature upwards, a 
discretization effect which here, since the lattice spacing 
is just the inverse lattice extension, is also a finite size 
effect. Furthermore, for small ll one expects to see finite 
temperature corrections according to l|5|). 

In order to study these finite temperature corrections 
systematically, we analyzed lattices with short tempo- 
ral extent, i.e., lattices with non-vanishing temperature. 
Fig. 2 shows the corresponding results, where we again 
plot the subtracted energy density as a function of fx . 

The lattice with the shortest temporal extent, 128 3 x 8, 
which corresponds to the largest temperature, shows a 
clear curvature. This curvature is due to the T 2 ll 2 /2 
term in (O, which appears as a square root when plotted 
as function of ll 4 . The effect is visible also for the other 
lattices, but becomes less pronounced as the temporal 
extent is increased, i.e., the temperature T is lowered. 
In order to study this effect quantitatively, we fit the 
finite temperature results to the continuum form §5§ plus 
two terms even in ll which parameterize the cutoff effects 
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FIG. 2: The energy density e(/i) — e(0) as a function of /i 4 , 
now for finite temperature lattices (all in lattice units). 

observed in Fig. 1. The fit function is given by 

C 2 Ll 2 + C A LL A + C 6 Ll 6 + C$Ll 8 . (14) 

Due to ([5]) the coefficient of the quadratic term should 
scale with the temperature such that one expects 

c 2 ~ T 2 /2 = 7V 4 " 2 /2 . (15) 

The coefficient for the quartic term should be constant, 

c 4 ~ 1/4tt 2 = 0.02533. (16) 

The results of the fit for the data used in Fig. 2, and for 
the largest lattice of Fig. I are given in Table 1. 

The table shows that with increasing N4 the two phys- 
ically significant parameters C2 and C4 approach the val- 
ues expected from the continuum formula ([5]): ci gets 
closer to N^ 2 /2 as listed in the second column, and C4 
approaches I /47r 2 = 0.02533. For the largest finite tem- 
perature lattice 128 3 x 24 the discrepancy is down to 9 % 
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C2 


C4 


C6 


C8 
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0.007812 


0.010125 


0.03519 


0.010 


-0.021 


12 


0.003472 


0.004125 


0.03178 


0.023 


-0.013 


16 


0.001953 


0.002192 


0.02803 


0.029 


-0.015 


24 


0.000868 


0.000947 


0.02587 


0.025 


-0.030 


128 


0.000030 


0.000032 


0.02543 


0.015 


0.016 



TABLE I: Results of the fits to the form JT4J). The spatial 
volume is always 128 3 . The temporal extension JV4 is given 
in the first column. In the second column we list the corre- 
sponding value of N^ 2 /2 which is what one expects for the 
fitting coefficient C2 in the third column. The coefficient C4 in 
the fourth column is expected to approach the constant value 
1/4tt 2 = 0.02533. 
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FIG. 3: The ratio — e(0))/ii 4 as a function of fj, (in 

lattice units). We compare the results for overlap to those 
from Wilson fermions. 

for C2, and 2 % for C4. The larger discrepancy for small 
N4 can be understood as a discretization effect, since the 
temporal lattice spacing 04 is related to the temporal ex- 
tension through 04 = I/-/V4 and thus larger N4 implies a 
smaller 04 . For comparison we also display the fit results 
for the 128 4 lattice, which corresponds to zero temper- 
ature. There we find excellent agreement (less than 1% 
discrepancy) for the parameter C4, governing the leading 
term at T = 0. The overall picture obtained from the fit 
results is that overlap fermions with chemical potential 
reproduce very well both, the \x term, as well as the fi- 
nite temperature contribution T 2 /i 2 /2. We conclude that 
the analytic continuation of the sign function does not in- 
troduce lattice artifacts, such as the /J? /a 2 term known 
to be present in a naive implementation of the chemical 
potential. 

In the final step of our analysis we study the discretiza- 
tion effect for larger values of fi and compare the re- 



sults to the data from the standard Wilson operator. In 
Fig. 3 we plot the ratio (e(/i) — e(0))/[i 4 as a function of 
(i. In the continuum at T = this ratio has the value 
1/4.7T 2 = 0.02533 indicated by the horizontal line. For 
small fj,, up to about fi ~ 0.7, the Wilson and overlap data 
fall on top of each other. For very small /1 both opera- 
tors show a prominent increase which is a left-over finite 
temperature effect, which for the ratio (e(/i) — e(0))//u 4 
shows up as a l/^ 2 term. In the range between fi = 0.1 
and 0.5 the data are close to the continuum value. Be- 
yond 0.5 the discretization effects kick in and the overlap 
and Wilson results start to differ. A comparison with 
the equivalent plot in where the results from various 
other lattice Dirac operators were presented, shows that 
the discretization effects of the overlap operator at large 
/j, are comparable to other formulations. 

V. SUMMARY 

In this article we have analyzed the energy density of 
the overlap operator at finite chemical potential. Follow- 
ing [i[ , the sign function in the overlap was implemented 
through the spectral theorem using the analytic continu- 
ation of the sign into the complex plane. The subtracted 
energy density e([x) — e(0) was analyzed for finite and 
zero temperature lattices. Fits of the data show that the 
expected continuum behavior is approached. No trace 
of unphysical n 2 /a 2 terms was found. We conclude that 
overlap fermions with chemical potential [i[ provide both, 
chiral symmetry and the correct description of fermions 
at finite density. 
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